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We use the internal- variable, effective-temperature thermodynamics developed in two preceding 
papers to reformulate the shear-transformation-zone (STZ) theory of amorphous plasticity. As 
required by the preceding analysis, we make explicit approximations for the energy and entropy 
of the STZ internal degrees of freedom. We then show that the second law of thermodynamics 
constrains the STZ transition rates to have an Eyring form as a function of the effective temperature. 
Finally, we derive an equation of motion for the effective temperature for the case of STZ dynamics. 
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I. INTRODUCTION 

Understanding the irreversible deformation of amor- 
phous systems remains a major challenge in nonequi- 
hbrium statistical physics and materials science P, Q. 
Systems of interest include noncrystalline solids below 
or near their glass temperatures, dense granular mate- 
rials, and various kinds of soft materials such as foams, 
colloids, and the like. An ongoing effort to develop a dy- 
namical theory of such systems has been based on the 
shear-transformation-zone (STZ) model of 0. Recent 
work has extended the original model to include an ef- 
fective disorder temperature as an essential ingredient 

BaaaiisQ. 

Our main goal in this paper is to develop an STZ the- 
ory that is consistent with the internal- variable, effective- 
temperature thermodynamics described in two preceding 
papers [H, [T^ . In [ill we focused on the role of inter- 
nal state variables in determining the nonequilibrium dy- 
namics of amorphous, not necessarily glassy, systems. We 
used the statistical interpretation of the first and second 
laws of thermodynamics to obtain equations of motion 
for the internal variables, and we emphasized the need 
to understand how both energy and entropy are shared 
between the internal variables and other degrees of free- 
dom. 

In jl^ we extended this development to include an ef- 
fective disorder temperature. Our basic premise in that 
paper was that the slow configurational degrees of free- 
dom of such materials are only weakly coupled to the 
fast kinetic/ vibrational degrees of freedom, and there- 
fore that these two subsystems can be described by dif- 
ferent temperatures during deformation. Using the tools 
of nonequilibrium statistical thermodynamics, we derived 
a general form for the equation of motion for the effec- 
tive temperature, and obtained a set of second-law con- 
straints on the thermomechanical equations of motion for 
such systems. 

We start here in Sec. HI] by summarizing the major re- 
sults of [l^, [l^ in a form appropriate for the STZ analy- 
sis. In Sec. mil we introduce the STZ degrees of freedom 
as thermodynamically well defined internal state vari- 



ables with associated energies and entropies. We then 
deduce specific forms for the STZ equations of motion 
based on the thermodynamic analysis. Our most impor- 
tant departure from earlier versions of the theory is that 
the STZ transition rates are now required to have an 
Eyring form as a function of the effective temperature 
rather than the reservoir temperature. In Sec. IIV[ we 
discuss the noise strength that determines the STZ anni- 
hilation and creation rates, and we derive an equation of 
motion for the effective temperature. Section IVl contains 
a summary of the STZ equations. We conclude in Sec. 
IVII with remarks about the significance and limitations 
of this theory. 



II. THERMODYNAMIC CONSTRAINTS 

We consider the deformation of an amorphous material 
in contact with a thermal reservoir at temperature On. 
We assume that Ofi is either below or not too far above 
the glass temperature 6g, so that the two-temperature 
theory developed in [3l is applicable. We express tem- 
peratures in units of energy, and set Boltzmann's con- 
stant Ub equal to unity. For simplicity, we assume from 
the beginning that the system is spatially uniform and 
that it undergoes only volume-conserving, pure-shear de- 
formations. 

The total, extensive, internal energy of this system, 
including a thermal reservoir, is 



Uc{Sc,E,i, {A„}) + Uk{Sk) + UrISr) , (2.1) 



where Uc and Uk, respectively, are the configurational 
and kinetic/ vibrational internal energies, Sc and Sk are 
the respective entropies, and {Aq} denotes a set of inter- 
nal state variables, soon to be identified as the STZ vari- 
ables. Ur is the energy of the thermal reservoir, which we 
assume to be strongly coupled to the kinetic/ vibrational 
subsystem. Egi is a deviatoric (traceless, symmetric), 
elastic shear strain. Note that our assumption of volume- 
conserving, pure-shear deformation allows us to omit any 
volume dependence in Uk, cf. Eq. (3.1) in [l3| . 
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The effective temperature x and the ki- 
netic/vibrational temperature 9 are 



X ■ 



dUc 
dSc 



fdUK 



K/v, 



(2.2) 



We assume that O^Oji, i.e. that the kinetic/ vibrational 
subsystem is always in equilibrium with the thermal 
reservoir. The shear stress acting on the configurational 
subsystem is 



Vsc = 



dUc 



Sc,{A»} 



(2.3) 



where V is the fixed total volume. As explained in Il3'|, 
the kinetic/ vibrational subsystem has no shear modulus, 
but it can support a viscous stress in the presence of 
shear flow. For further simplicity, we assume that the 
kinetic/ vibrational viscosity vanishes. 
The total entropy is 

Stot ^ SciUc, -Ee,, {A,}) + Sk{Uk) + Sr{Ur) . (2.4) 

The expression for any one of these three entropies can 
be inverted to obtain the corresponding internal energy 
function in Eq. (j2.ip . or vice versa. 

Without further loss of generality, we specialize to the 
case of pure, planar shear oriented along fixed axes, say x 
and y, and define sc = sc,xx = — sc,yy We assume (for 
small elastic deformations) that the rate of deformation 
tensor is the sum of elastic and inelastic parts, D = Dei + 
Dm, where Dg; = Ee,i, and we define An = Din,xx = 



-A. 



All other elements of these deviatoric tensors 



vanish; thus, for example, the rate of inelastic work done 
by the shear stress is sc ■ D.^ = 2sc Din ■ 

The analysis in Ts] produced an equation of motion for 
the effective temperature that is basically a statement of 
the first law of thermodynamics, i.e. a heat-flow equa- 
tion. For the present case, this equation has the form 



a 



Here 



eff 



X = >Vc(sc,{A4) + A(x,0) . (2.5) 



Cf X = xSc 



(2.6) 



is the time rate of change of the heat of configurational 
disorder, and Cy^ is an effective (extensive) heat capacity 
at constant volume. As in the preceding papers, the non- 
negative dissipation rate Wc - the difference between 
the rate at which inelastic work is being done on the 
configurational subsystem and the rate at which energy 
is being stored in the internal degrees of freedom - is 



Wcisc, {Aa}) ^2Vsc D,n-J2 



fdUc 



Sc,E,. 



Aa > 0. 

(2.7) 

Non-negativity of Wc is an important second-law con- 
straint that plays a central role in the analysis to follow. 



The second term on the right-hand side of Eq. (|2.5p , 

A(x,0)(l-|)^Q (2.8) 

is the rate at which heat is fiowing into the configura- 
tional subsystem. Here, A{x, 0) is a non-negative ther- 
mal transport coefficient that, as will be seen, depends 
on other dynamical variables in addition to x and 0. 

III. STZ EQUATIONS OF MOTION 

The basic assumptions of the STZ theory have been 
described in To the extent possible, the following 
discussion follows the steps outlined in that paper. 

The main idea is that deformation of amorphous mate- 
rials occurs via localized molecular rearrangements that 
take place at shear transformation zones (STZ's). The 
STZ's are created and annihilated either by thermal fluc- 
tuations or by noise generated by the deformation it- 
self. They are rare, ephemeral fluctuations that are es- 
pecially important for irreversible deformations because 
they make stress-driven transitions between two, energet- 
ically almost degenerate orientations. Thus, the STZ's 
are two-state systems. There is nothing arbitrary about 
this two-state picture. The STZ's have the special prop- 
erty of being able to shift between one orientation and 
another in response to a shear stress. Sites with this 
property are already statistically unlikely, and higher- 
order degeneracies are statistically negligible. 

The difference between what we are doing here and 
the analysis presented in 0] is that now, on the basis of 
[13, [l^, we insist on a proper thermodynamic descrip- 
tion of the STZ's as internal degrees of freedom. Such a 
description requires a specific STZ model. To construct 
any such model, we must make physical assumptions that 
may need to be modified in later applications. In partic- 
ular, as in the earlier work, we assume that there is just a 
single kind of STZ, with a single characteristic formation 
energy ez, and a single mechanism for making transitions 
between the two orientational states. 

For additional simplicity, we go back to the original 
version of the theory 3] in which the STZ's occur only 
with orientations either "+" or "— " with respect to the 
shear direction. A procedure for averaging over STZ ori- 
entations and constructing a properly invariant tensorial 
version of the theory was presented in That proce- 
dure works just as well for the present analysis, but seems 
unnecessarily complex for present purposes. The inter- 
nal state variables are the extensive numbers of STZ's in 
these two different states, and N^. As usual, define 



A = 



N 



N+ -N^ 
N+ + N^' 



(3.1) 



Thus, the set of internal state variables {Aq} reduces to 

Our arguments in |12l . [l3[ tell us that we must include 
the entropy associated with the internal variables A and 
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m in this analysis. If we take the two-state model liter- 
ally, then we compute this entropy by counting the num- 
ber of ways in which we can distribute N-^ "+" zones 
and iV_ "— " zones among, say, N available sites in the 
system. This number is 



exp {Sz) 



Nl 



(3.2) 



N+\N-liN - N+ - N-)] 
which, after use of Stirling's approximation, reduces to 
Sz{A,m) 



N 
where 



So{m) = In 2 (1 -f-?7i) ln(l + m) (1 — m) ln(l — m). 

^ ^ (3.4) 



To use this formula, write 

Sc^Sz{A,m) + SiiUi), 



(3.5) 



where Si and Ui, respectively, are the entropy and en- 
ergy of all the degrees of freedom of the configurational 
subsystem apart from those attributable to the STZ's. 
Accordingly, 

Uc{Sc,E,i,A,m) ^ N Aez + Ui{SuE,i) 

^NAez + Ui[Sc~SziA,m),E,i], (3.6) 

where ez is the formation energy of an STZ. Equations 
p.Sp and (|3.6p are equivalent to each other if we write 
Ui = Uc -NAez in Eq. (^3)1 . 

In terms of these STZ variables, the inequality in Eq. 
(12. 7p becomes 



^2Vsc D„,-N 

+NxA^m>0. 
am 



ez + X In 



A 



1 - A 



- xSoim) 



A 



(3.7) 



To make further progress, go back to the original STZ 
equations of motion for the N± 



To N± = Ri±sc) - RiTsc) N±+T [l-mi -N± 



(3.8) 

Here, tq is a time scale, the factors R{±sc)/tq are the 
rates at which STZ's switch back and forth between their 
two orientations, F/to is the rate factor for creation and 
annihilation of STZ's, and N'^'^ is an as-yet undetermined 
"equilibrium" value for the number of STZ's. The super- 
script "eq" is used here and below to denote steady-state 
equilibrium. Note that, in Eq. (|3.8p . we are assuming 
that the STZ creation rate is the same for both STZ ori- 
entations, independent of the orientational state of the 
system as a whole. 



The deviatoric, inelastic rate of deformation tensor is 
D^n = ^ [i?(+sc) - R{^sc) N+\ , 
where vq is a molecular-scale volume. As usual, define 
C{sc) ^ \ [R{+sc)+R{-sc)\, 



rise) 



Rj+sc) - Rj-sc) 
R{+sc) + R{-sc) 



-A InA-(l-A) ln(l-A)+A5'o(TO), (3.3) Then, 



To Din = ^^^^ AC(sc) [T{sc) - m] . 



(3.9) 



(3.10) 



In previous papers, we defined N vq/V = Cq. We will re- 
turn to this notation in Sec. |Vl The equations of motion 
for A and m are 



roA = f(A^9- A), 
where A^« = N^i/N, and 

To m = 2C{sc) [r{sc) — to] — F to 



A 
A 



TO. 



(3.11) 



(3.12) 



The next step in this analysis is to impose the second- 
law constraint expressed in Eq. p.7p . We immediately 
encounter a difference between the present situation and 
the one described, for example, by Maugin in fT3|. Specif- 
ically, the inelastic rate of deformation I?i„ appearing in 
Wc is not simply proportional to the time derivatives 
A and to. Therefore, we cannot satisfy the inequality 
in Eq. (|3.7p by identifying the coefficients of those time 
derivatives as thermodynamic forces associated with en- 
ergy landscapes, and then requiring that A and to both 
relax toward free-energy minima. In fact, our situation 
is more interesting. It is almost certainly typical of open 
systems in which external work is being done and energy 
is being dissipated, and where no variational formulation 
is relevant. 

Our strategy is to use Eq. p.l2p to evaluate to in 
Eq. (|3.7p . and thereby to write Wc as the sum of two 
terms, one proportional to A, and the other proportional 
to the stress-dependent quantity T(sc) — to. These two 
terms must individually be non-negative. The inequality 
in (|3.7p becomes 



^Wc(^c,A,TO) = -fxATO^ 



X In 



A 



1 - A 



X Ao(to) +Xm — 



(3.13) 
To A 



-2AC{sc)[T{sc)-m] (^vosc + x^^ 



> 



From Eq. (|3.4p . we know that 

dSa 1 , / 1 + to\ , -1 / n ,x 

= -- In = -tanh ^(to). (3.14) 

dm 2 V 1 — TO / 
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Therefore, the first term in the expression for Wc in Eq. 
p.l3p is always non-negative, and we can set it aside for 
the moment. 

The second term in Eq. (|3.13p produces a standard, 
variational, second-law inequality of the form 



dFz 
dA 



A > 0, 



where 



Fz{A,m)=NezA-x 



SziA, m) — m 



dSz 



dm 



(3.15) 



(3.16) 



is a free energy. A*^* in Eq. (|3.1ip must be the value of 
A at which 



Therefore, 



where 



dFz 
dA 



A'\x,m) 



= 0. 



A=A"i 



2;eq 



Z'^'^{x, m) = exp 



h JQ(m} — m — — 

X dm 



(3.17) 



(3.18) 



(3.19) 



For X ez, we expect A*^"? w Z'^^ <C 1, which is con- 
sistent with the basic idea of a low density of STZ's. 
We then obtain the expected Boltzmann factor, A"^"* k, 
exp(— e^/x), with a small modification from the m- 
dependent entropy. The term proportional to m dSg /dm 
in Eq. (|3.19p means that Z*^"^ diverges weakly, and 
/^eq _^ 2^ when TO ^ ±1. However, it is easy to see from 
the denominator in the equation of motion for m, i.e. ei- 
ther Eq. (|4.6p or Eq. (|4.8p shown below, that to — > ±1 is 
a dynamically inaccessible limit. Therefore, so long as 
ez is the largest energy scale in the problem - which has 
always been the case in prior applications - the require- 
ment of small A is satisfied. 

The more interesting result comes from the term pro- 
portional to T{sc) — TO in Eq. ((3T3)) . That term must 
be non- negative for all values of the stress sc, i-e. 



[r{sc) 



«osc + X^)>0, (3.20) 
dm 



which means that the two stress-dependent factors, 
T{sc) — m and vq sc + xdSo/dm, must each be monoton- 
ically increasing functions of sc that change sign at the 
same point for arbitrary values of to. From Eq. (|3.14p . 
we see that this condition can be satisfied only if 



T(sc) — tanh 



^■0 Sc 
X 



which, according to Eq. (|3.9p . means that 



X 



Rise) = RoiscX,^) exp 



(3.21) 



(3.22) 



where Rq is a symmetric, non-negative function of sc- As 
indicated, Rq may also depend on the temperatures x and 
9, because the transitions between STZ orientations are 
very likely to be thermally activated processes. Equation 
([3^ indicates a major difference between the present 
thermodynamic results and the earlier theories. In the 
latter, we started with physical models for the transition 
rates i?(±sc'), and then assumed that the dependence of 
the internal energy on the STZ variables would be con- 
sistent with these rates. Here we start with a known 
internal energy, and must argue in the other direction to 
make sure that the rates are consistent with thermody- 
namics. In particular, Eq. ((3?22)) tells us that the STZ 
transition rates must have an Eyring form with the effec- 
tive temperature x rather than the reservoir temperature 
9 in the exponent. 



IV. NOISE STRENGTH AND EQUATION OF 
MOTION FOR X 

Having used the second law to deduce equations of 
motion for the STZ variables, our next steps are to go 
back to the first law in Eq. (|2.5p and use the expressions 
for A and to to compute F, and then to derive the STZ 
version of an equation of motion for x- Both of these 
steps again require going beyond purely thermodynamic 
arguments, and making additional physical assumptions. 

Equation (|2.5p now can be expressed explicitly in terms 
of the internal variables: 



cfx~A{x,e) (1 



X 



xSc 



(4.1) 



-N 



To 

2NA 



tq oA V ^ To dm 

dSo 



To 



C{sc) [r{sc) 



vqSc + X 



dm 



As in previous STZ papers, we assume that the rate 
factor F is a sum of two independent noise strengths, 
F = T{sc, x) + p{(^)- Here F(sc, x) is the part of the rate 
factor determined by mechanically generated noise, and 
p{9) is the super-Arrhenius, thermally generated part. 
We next invoke Pechenik's hypothesis [15|, which iden- 
tifies F as being proportional to the total rate of heat 
production per STZ 



xSc-Q = Wc ^ ^('c,x) jvA^io^o, 



To 



(4.2) 



where the proportionality factor sq has the dimensions of 
stress. Inserting this relation into Eq. (|4.ip and solving 
for F, we find 



Tisc,x)+P{0) 



A/"(sc, A,m) 
A(A,to) ' 



(4.3) 



where 



A/'(sc, A, to) = p{9) vo So 
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+2 C(sc) [T{sc) - H (^0 sc + X ^) (4.4) 



and 



A(A,H = ^o.o + ^ — -1 Urnx^. (4.5) 



The equation of motion for m, Eq. (|3.12p . becomes 

>1(sc, A,to) 



Torn 



A(A,m) ' 



(4.6) 



where 



A4(sc, A, m) = - TO — ^(6*) sq I'd 



+ 2C(sc) [T(sc)-™] 



Wo (so ~ msc) 



"a 



- 1 



dSn OF. 



dm dA 



(4.7) 



At this point, it is useful to distinguish between slow 
and fast processes, as was done in [a, Q- The inelastic 
deformation rate given in Eq. (j3.10|) contains a factor A, 
meaning that it is proportional to the density of STZ's 
and is small. The equation of motion for x will be seen 
to be similarly slow. On the other hand, the equations of 
motion for A and to contain no such factors A. These in- 
ternal state variables respond rapidly to changes in their 
environments. Therefore, we simplify the analysis by set- 
ting A = A*^"?, and replacing m by m^'^, the stationary 
solution of 



To TO 



2C{sc) [T{sc) - to] (l - m sc/so) - m pi^) 



I + (mx/voso) {dSo/dm) 

(4.8) 

This solution is shown explicitly in Eq. (|5.5p . These ap- 
proximations are always valid for steady-state solutions 
but, as seen in they also work well for transients. 

In steady state, and at low temperatures where p{9) « 
0, Eq. (jTS]) exhibits the usual [1, S Ell exchange of 
stability at a yield stress (minimum flow stress) Sy deter- 
mined implicitly by 



Sy tanh 



^'0 Sy 

Xa 



So, 



(4.9) 



where xo is the steady-state value of x in ti^c limit of 
vanishingly small strain rate. According to Eqs. (|4.8p 
and (jS.Sp . for p{9) = 0, m*^^ goes through a maximum 
value of tanh (vq Sy/xo) sc = Sy. At that point, Eqs. 
(PrS)) and (|XTg)l tell us that the condition 

A^«(xo, m^«) « exp (~ ££^^) « i (4.10) 
V Xo J 

requires that ez be much larger than xo and vq Sj,, which, 
as noted earlier, is generally true. 



To complete this development, we need an explicit 
equation of motion for x, and again we need to make ad- 
ditional physical assumptions. Use Eqs. (|4.ip and (|4.2p 
to write 

Cfx^^^^^^NA^^voso + Aix,e) (l-f). 

(4.11) 

The thermal transport coefficient A{x, 0) is one of two 
places in this theory where the weak coupling between 
the configurational and kinetic/ vibrational subsystems 
must be modeled explicitly. The other place is the noise 
strength F defined in Eq. (|4.2p . where we argued that me- 
chanically generated noise contributes additivcly, along 
with the thermal noise, in creating configurational disor- 
der. Similarly, it seems plausible that the overall heat 
exchange between the two subsystems is enhanced by 
mechanical noise. Thus we propose that A have a form 
similar to that of F, and write 

A(x,0)-^^[F(sc,x) + ^p(^)], (4.12) 

To 

where k is a dimensionless parameter, the factor 9 has 
been inserted for dimensional reasons, and ao is a dimen- 
sionless quantity to be determined as follows. 

Separate the right-hand side of Eq. ()4.1ip into parts 
proportional to F and p, and then write this equation in 
the form 



To 



X 



N 



= T{sc,x) [A"" vo So + ao {9 - x) 
+ aoKp{9){9^x)- 



(4.13) 

In it was argued that athermal (p — 0) amorphous 
systems reach steady state for effective temperatures x 
equal to some function x{q), where g is a dimensionless, 
non-negative measure of the total strain rate. For time- 
independent stresses, q is the magnitude of To Din- This 
means that the quantity in square brackets in Eq. (j4.13p 
must vanish at x = x{q), a condition that we satisfy by 
setting 



ao 



A'l VQ so 



(4.14) 



Thus, Eq. ((4T3)) becomes 



vo So N 
A^i 

x{q) - 



(4.15) 



[r(sc,x) im-x)+^ip{0) {9-x)] 



Equation (|4.15p is essentially the same x equation that 
we have used in previous applications. The main differ- 
ence is the prefactor (x — 9)^^. Non- negativity of ao 
requires that xil) > which is a plausible and interest- 
ing constraint. The steady-state solution of Eq. ()4.15p 
is 



X^ 



T{sc)xiq)+p{0)t 
risc)+pi9) 



(4.16) 
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The function r(sc) vanishes in the limit of vanishing 
strain rate q; therefore, for fixed, nonzero p{0), Xss 
as g ^ 0. On the other hand, if the strain rate is fixed 
and p{6) becomes small, then Xss x{q)- As pointed 
out in , the crossover between these limiting behaviors 
takes place at very small strain rates for small p{0), and 
therefore it can be very difficult to determine whether a 
glass transition has occurred. At higher temperatures, 
this crossover occurs at higher strain rates, and the con- 
dition x{q) > 9 requires that x be a function of in some 
circumstances. For the moment, we note that physically 
realistic systems do not probe the extreme limit of van- 
ishingly small strain rate, and we therefore assume that 
x{q) — 9 ^ Xo ^ is a, positive constant for situations in 
which the system is deforming at experimentally accessi- 
ble rates. 



V. SUMMARY OF STZ EQUATIONS 

We conclude this part of the paper by summarizing 
the STZ equations in their most usable versions, that is, 
in the limit in which the relaxation of the STZ variables 
A and m is much faster than the rates at which plas- 
tic deformation and the effective temperature respond to 
changes in the external driving forces. Many of these 
equations are the same as the ones that appear - in more 
general tensorial versions - in As noted previously, 
however, there are some differences. 

The rate of inelastic deformation, given here in Eq. 
(I3.10p . is a function of the configurational shear stress 
sc (assuming no appreciable contribution from the vis- 
cous stress in the kinetic/vibrational subsystem) and the 
effective temperature x 



and 



where 



-ez/x 



(5.1) 



(5.2) 



T = tanh 



X 



(5.4) 



Ro{sc) is an arbitrary, symmetric function of the shear 
stress Sc- rn'^''{sc,9) is the stationary solution of Eq. 



m^'(sc,e) 



50 

2sc 



f£r- 

So 



2C 



So 
'2sc 



sc_ 
So 



2C 



4^r. 

So 



(5.5) 



The parameter sq is a stress that can be determined from 
the low-temperature yield stress (minimum flow stress) 
Sy via Eq. ([49]) 



Sy tanh 



vq Sy 

Xo 



= So, 



(5.6) 



where xo is the steady-state value of x in the limit of 
vanishingly small strain rate. 

It is useful to look at the equation of motion for x, Eq. 
(|4.15p . in two special cases. First, consider the param- 
eter range relevant for deformations of ordinary plastic 
materials such as metallic glasses. The experience gained 
from the studies reported in [7] and [1] suggests, for tem- 
peratures not too far above the glass transition, and for 
strain rates not extremely small, that we can assume that 
x{q) « Xo remains constant at a value larger than 9, so 
that the dimensionless quantity x/(x — 9) is a slowly 
varying function of 9 that can be absorbed into other pa- 
rameters such as the effective heat capacity and k. When 
this is true, Eq. (|4.15p can be written in the form 



ToCoX = e '^^/'^ 



Tiscx) [l-^^)+kp{9) (l-l 



, (5.7) 



and 

/(sc,x) = -C(r- m'^'^). 

To 

Here, we have reverted to the earlier notation, eo = 
N vq/V, which is the ratio of a molecular volume vq asso- 
ciated with STZ transitions to the volume per molecule 
in the system as a whole, and is of the order of unity. 
The STZ formation energy ez previously was denoted 
hy ksTz- In [7], Tz was found to be larger than the 
glass temperature by a factor of about 30 for a metallic 
glass; and the time constant tq was of the order of a fem- 
tosecond. We have abbreviated the functions C and T as 
follows 



C = Ro{sc,X,(^) cosh 



^■0 Sc 

X 



(5.3) 



where cq and k are dimensionless constants of the or- 
der of unity. To use this equation, we need the explicit 
expression for F 



r(sc,x) 



where 



^/■(scx) 



(m^^x/so^o) tanh \m<"i)' 



(5.8) 



A^(sc,x) = p{0) -t&n\i-\m^'^) 

vo So 

+2C(sc) [T{sc)-m''i] 



Sc X 



So Vq So 

and m'^''(sc,x) is given by Eq. 



tanh-\m'=«) 



(5.9) 
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Second, consider the athermal limit of Eq. (|4.15p by 
setting 9 = and p{9) = 0. In this case, we have 

TO So X = e~ ^-/^ r(sc, x) (^1 - ^) , (5.10) 

where, now, cq = Cy-^ /eo, and Cy^ is the effective heat 
capacity per unit volume in units of Boltzmann's con- 
stant kB- This limit is appropriate for granular materi- 
als, bubble rafts, and the like, where ordinary thermal 
fluctuations are irrelevant, and the disorder described by 
the eff'ective temperature is generated only by externally 
driven deformation. Thus, only states with stresses above 
the yield stress are relevant, and Eq. (|5.5p tells us that 
m'^'^ = Sq/sc (exactly). Moreover, when sc ^ sq, we 
have 

soA^9r(sc.,x) «2scA«, (5.11) 

so that the noise strength is just proportional to the rate 
at which inelastic work is done on the system. We have 
used x{q) on the right-hand side of Eq. (|5.10p . instead 
of its small-g limit xo, because large values of q are more 
easily attainable for systems in which the intrinsic relax- 
ation time To is not microscopically small. As shown in 
[TtI . x{q) increases rapidly when q grows to values of the 
order of unity. Thus the restoring term in Eq. (|5.10p 
becomes small; and the resulting rapid growth of x pro- 
duces localized shear failure. This mechanism was shown 
in [isj to provide a plausible explanation of rapid stress 
drops and localized failure in earthquake faults. 



VI. CONCLUDING REMARKS 

We have made many simplifying assumptions in devel- 
oping this thermomechanical version of the STZ theory. 
Some of these assumptions were needed only to simplify 
the presentation, and seem to have little if any physi- 
cal importance. For example, it should not be difficult to 
rewrite this theory in tensor notation, as in f^, and apply 
it to spatially nonuniform situations with orientationally 
varying stress and flow fields. It will be technically more 
difficult to deal with situations in which both volumetric 
and shear deformations are occurring and are coupled to 
each other; but here again there seems to us to be no 
problem in principle. 

Yet another example of simplification is that, through- 
out this series of three papers, we have dropped terms 
that would have described thermoelasticity or, more per- 
tinently in the context of nonequilibrium phenomena, 
thermo-viscoelasticity. Here too, we see no intrinsic dif- 
ficulties. In fact, we see attractive opportunities to use 
a thermo-viscoelastic version of this theory for studying 
the behavior of glasses subject to thermal cycling in the 
neighborhood of the glass temperature. 

One of our more problematic simplifications is our as- 
sumption that we can distinguish elastic from plastic 



strains, and use the elastic strain as an independent ar- 
gument of thermodynamic functions such as the internal 
energy or the entropy. As we have stated here and in ear- 
lier papers, we maintain that the plastic strain, necessar- 
ily measured from some reference configuration (possibly 
evolving) , cannot be a physically meaningful variable for 
determining the current state of the system or predict- 
ing its subsequent motion. Thus, we have insisted on 
expressing our equations of motion in Eulerian coordi- 
nates, and using the internal state variables to carry the 
memory of recent deformations. 

This self-imposed requirement leaves us with an as- 
yet unsolved problem regarding elasticity. The problem 
is compounded here by our recognition of the extended 
thermodynamic roles played by internal degrees of free- 
dom, which, as we have seen, may store energy in recov- 
erable forms as well as relax irreversibly toward states 
of equilibrium. In such situations, it is unclear to us 
whether "elastic" behavior is always the same as "re- 
versible" behavior, or whether the conventional Kroner- 
Lee (lol . [20j decomposition of elastic and plastic displace- 
ments is generally correct. We have evaded these issues 
so far by restricting our attention to infinitesimally small 
elastic displacements. However, we suspect that these 
questions now require more serious attention. 

Our list of topics needing further investigation includes 
the choice of rate factors in the STZ theory. Our most no- 
table departure from earlier STZ results is the relatively 
simple, x-dependent transition rate shown in Eq. (j3.22p . 
This formula is primarily a result of our statistical inter- 
pretation of the second law of thermodynamics in p^ : 
it is related to the two-temperature theory only in the 
sense that it is the effective temperature x, and not the 
thermal temperature 0, that governs the configurational 
subsystem's motion toward statistically more probable 
states. So far as we can tell, this result does not sub- 
stantially change previous conclusions, e.g. in 0, In 
fact, the stronger stress dependence in Eq. p.22p may 
be needed in order to understand seismic data [21]. 

This statistical interpretation of the rate factors is 
especially difficult for jammed states at low tempera- 
tures, where the stress is below the yield stress and 
p{0) = 0. Our theory predicts that, in this situation, 
m = tanh(wo sc/x)- This result makes sense for a glass 
below its glass transition temperature, where thermal 
fiuctuations still can activate transitions between the 
states of STZ's even if they cannot create new ones. In 
this case, we can change the inelastic strain by chang- 
ing the stress, although reequilibration to a new state of 
deformation might be very slow. 

For a granular material, however, the most we can 
say is that m = tanh(uoSc/x) is the statistically most 
likely average orientation of STZ's at the given values 
of Sc and x- Such a state might be achieved by tap- 
ping the system, i.e. by artificially introducing some- 
thing like thermal noise. But the way in which such a 
jammed system responds to changing stresses has to do 
with whether it forms force chains or bridging structures 
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or the like. Such mechanisms cannot be included in a 
theory of the kind we are discussing here. Therefore, 
when talking about granular materials in Sec|Vl we have 
restricted ourselves to unjammed systems that are un- 
dergoing deformation. More generally, this limitation of 
the STZ theory emphasizes the need for a more thorough 
investigation of the limits of validity of this theory and of 
similarly constructed statistical theories of noncrystalline 
deformation. 
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